The Mesoscopic SNS Transistor 
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In a mesoscopic superconductor - normal metal - superconductor (SNS) heterostructure the quasi- 
particle distribution can be driven far from equilibrium by a voltage applied across the normal metal. 
This reduces the supercurrent between the superconducting electrodes, which creates the possibility 
to use these SNS junctions as fast switches and transistors. We describe the system in the frame- 
work of the quasiclassical theory and find good agreement with recent experiments. We propose 
further experimental tests, for instance the voltage-dependence of the current-phase-relation, which 
includes a transition to a 7r-junction. 



the normal wires between the normal reservoirs is shorter 
than the inelastic scattering length, the quasiparticle dis- 
tribution function / has the noncquilibrium (twa,step) 
form observed in the experiments of Pothicr ct ala. 
(iii) The noncquilibrium distribution influences (reduces) 
the supercurrent between the two superconducting con- 
tacts through the normal metal. 



Nonequilibrium effects in superconducting systems 
have been gaining new attention due to the increased 
activities in the field of mesoscopic electron transports. 
In contrast to earlier work on nonequilibrium supercon- 
ductivity the new experiments show non-local and size- 
dependent (d) effects by reaching temperatures below the 
characteristic Thouless energy E^h = T>/d 2 . The under- 
standing of this regime is not only of fundamental inter- 
est, but also important for nano-electronic applications. 
The size reduction of electronic devices is accompanied 
by an increase in operation frequency. For instance, in 
the system considered below the latter is limited by i?Th- 

Recently Pothier et ala probed the quasiparticle prop- 
erties in short diffusive wires by coupling tunneling con- 
tacts to it. They found that in mesoscopic wires the 
distribution function has a nonequilibrium energy depen- 
dence, with a double-step structure at the electrochem- 
ical potentials of both reservoirs. From the smearing of 
this distribution they further drew conclusions about in- 
elastic relaxation processes. In different setups, sketched 
in Figs. 1, Morpurgo et alB demonstrated that in a super- 
conductor - normal metal heterostructure the nonequi- 
librium quasiparticle distribution due to a normal current 
flow in N can be used to tune the supercurrent between 
the superconducting electrodes. This opens the perspec- 
tive to use such devices as ultrafast transistors. 

Topaccount for their experimental findings, Morpurgo 
et ala proposed a qualitative model based on a quasi- 
equilibrium distribution function with locally enhanced 
effective electron temperature. This picture is appropri- 
ate in the limit of strong electron-electron interactions. 
However, inelastic processes have only a weak effect in 
the mesoscopic sample considered here. In fact we find 
the best performance of the devices in the opposite limit. 

In this article we will describe mesoscopic supercon- 
ductor - normal metal heterostructures as shown in Fig. 
1 a and b. For a quantitative analysis we use the quasi- 
classical theory. It accounts well for the relevant physics: 

(i) The spectral properties in the normal metal are mod- 
ified by the proximity effect due to the presence of the 
superconducting electrodes. 

(ii) The nonequilibrium quasiparticle distribution func- 
tion is found as the solution of a kinetic equation. If 
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FIG. 1. Different realizations of the SNS transistors as used 
in Refs. [|B|]. The supercurrent is tuned by (a) a perpendicular 
or (b) a parallel normal current flow. 

For best performance as a transistor, in the setup of 
Fig. la, the width of the superconducting contacts ds 
should be chosen narrow compared to the width and 
length of the normal wire, d and L, respectively. The 
first condition assures that only a small fraction of the 
control normal current is diverted through the supercon- 
ductors, and accordingly the quasiparticle distribution 
function in N is little disturbed by the presence of the su- 
perconducting electrodes. The second condition assures 
that the voltage is nearly constant along the supercon- 
ducting leads, while the total voltage drop responsible 
for the nonequilibrium and reduction of the supercurrent 
may be large. Since the effect relies on a deviation of the 
distribution function from local equilibrium with shifted 
electro-chemical potential, the normal wire should have 
mesoscopic dimensions, i.e. the length L should not ex- 
ceed the inelastic relaxation length l m . 

The presence of the superconducting electrodes in- 
duces correlations in the normal metal (proximity ef- 
fect), which are responsible for a supercurrent. Their 
decay length depends on energy: correlations with en- 
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ergy e 3> Etyi decay exponentially, while those within 
the range of order i?Th carry the supercurrent. If, in a 
nonequilibrium situation, these states are occupied and 
in this way blocked for superconducting correlations, the 
superconductivity is weakened and the supercurrent re- 
duced. This suppression mechanism will be described 
in the following, based on the real-time formalism of 
quasiciassical Green-Keldysh functions in the diffusive 
limitM. 

In the first step we describe the proximity effect in the 
normal metal by analyzing Usadel's equations. The stan- 
dard parameterization of normal and anomalous retarded 
Green functions G R = cosh a and F R = sinh ae lx allows 
us to write these equations, for the normal metal region 
between the superconducting electrodes, in the form 



Vdla = -2iesinha - (D/2) {d x \f sinh2a 
d x je = , j e = (d x x) sinh 2 a . 



(1) 



Here, T> is the diffusion coefficient and x the coordi- 
nate normal to the NS interfaces. We also introduced 
the energy-dependent 'spectral current' j e . For simplic- 
ity we ignored here a dependence of the spectral quanti- 
ties on the coordinate y parallel to the interfaces. This 
dependence should merely lead to a quantitative modi- 
fication of our results and conclusions. In the realistic 
limit A 3> i?Th and for transparent metallic interfaces 
the boundary conditions at these interfaces x = ±rf/2 
read (see [ 0] and Refs. therein) 



a(±d/2) = -in/2 , x(±d/2) = ±0/2. 



(2) 
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FIG. 2. The spectral current Im(j e ) as a function of energy 
for different values of the phase difference 4>. 

In two limits we find analytic solutions of Eqs. (Jj], ||). 
For low energies e <C 2?Th we obtain 

a ~ -iw/2 + (e/E Th ) a((f>) + 0(e 3 ) 
X ^0x/d-(e/E Th )H(0) + O(e 3 ) , 



where a(</>) and &(</>) are real-valued functions (omitted 
for brevity), while for high energies e 3> -Erh we obtairo 



F R =F (x- d/2)e^/ 2 + F (d/2 - x) e -^ /2 



(3) 



F (s) =M 



j(s) = i(V2- i) e -^V=2W^ 



, ... rl • , , (l + g 2 )(l + 6g 2 +g 4 ) , 
Im(j e ) = 64sm0Im ( (l-q^f m 



s=L/2 



In addition we have studied the problem numerically. 
Combining our results in Fig. 2 we observe the follow- 
ing features of the spectral current Im(j e ). It is an odd 
function of e and shows a proximity induced mini-gap£3, 
e g ~ 3.2i?Th at (j> — 0, below-which Im(j c ) = 0. This gap 
decreases with increasing c/O and vanishes at 4> — tt. At 
energies directly above the gap, Im(j e ) increases sharply, 
but rapidly decreases at higher e. At large energies, it 
changes sign and oscillates around zero with exponen- 
tially decaying amplitude. 

Next we determine the nonequilibrium quasiparticlc 
distribution in the normal metal between the reservoirs, 
which are at different electrochemical potentials ±eV/2. 
In a mesoscopic length wire L <C ^ n m the diffusive limit 
the distribution function obeys the kinetic equation 



dlf = 



(4) 



In the absence of superconducting contacts its solution 



(l/2-y/L)n(e + eV/2) 
+ {l/2 + y/L)r(e-eV/2) 



(5) 



has two temperature-rounded steps at the electrochemi- 
cal potentials of both reservoirs. The step heights depend 
on the position along the wire; in this way the distribu- 
tion function interpolates linearly between the boundary 
conditions at y = ±L/2. This functional dependence had 
been detected in the experiments of Pothier et ai.Q. 

Although the distribution function definitely has not 
a thermal form, a local electrochemical potential of the 
normal metal, fi(y), and an effective electron tempera- 
ture can be defined by the corresponding moments of 
the distribution function. For instance fi(y) follows from 
JZ.deif^y) - - - 0, where denotes 

the Fermi function. In the following we will consider 
the situation where the electrochemical potential of the 
superconductors coincides with the local value of the nor- 
mal metal, which guarantees that there is no net current 
out of the normal metal into the superconductors. Since 
we further have chosen the size of the superconducting 
contacts, ds, small compared to the width and length, d 
and L, of the normal wire, the distribution function in N 
is little disturbed by the contacts. 

On the other hand, the quasiparticle distribution func- 
tion influences the superconducting correlations induced 
in the normal metal. The supercurrent through the het- 
erostructure is given by 
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/s = 97r/ de[l-2f(e,y s )]lm(je), 



(0) 



where 1/Rd = 2e 2 No'DS/d, S is the junction area, and ys 
denotes the position of the superconducting electrodes. 
It is important to note that the energy is measured rel- 
ative to the electrochemical potential of the supercon- 
ductors, which we have chosen to coincide with the local 
value of the normal metal, fJ-(y). Due to the odd sym- 
metry of Im(j e ), only the odd component of a nonequi- 
librium quasiparticle distribution modifies 1$. Accord- 
ingly the first term in the integral (^) can be written 
asH —f(e,ys) + f(—e,ys), which displays that an excess 
number of electron-like or of hole-like excitations have 
the same effect on the supercurrent. 

The largest effect is found when the superconducting 
electrodes are placed symmetrically between the two nor- 
mal reservoirs (i.e. at y = 0, see eq. (^)). The resulting 
modification of the supercurrent across the SNS junction 
is presented in Fig. 3 as a function of the voltage V across 
the normal metal. 




FIG. 3. The supercurrent as function of control voltage 
T = 0) and temperature (at V = 0) for various values of q 
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At low temperatures T <§; eV, /(e, 0) deviates from the 
equilibrium value only in the window —eV/2 < e < eV/2. 
Since the spectral current Imj e vanishes for e < e g there 
is no modification of of the supercurrent for small volt- 
ages eV < 2e g . On the other hand, for eV > 2e g , extra 
quasiparticle and hole states with energies |e| below eV/2 
are occupied and the supercurrent is diminished. Since 
this energy window increases with increasing V the su- 
percurrent decays rapidly with voltage (cf . Fig. 3) . This 
is exactly what has been observed in the experimentsu. 
Furthermore, at still larger voltage eV > 10£/rh the su- 
percurrent changes sign since the integral in (|^) is dom- 
inated by the energy interval where hnj e < (F-ia;. 2). 
We thus find a transition to a so-called 7r-junctionE3, con- 
trolled by nonequilibrium effects. This effect is rather 
pronounced, the critical current of the 7r-junction is ap- 



proximately 30% of Z c at T = eV = 0. 

For high voltages or temperatures eV, T ^> e g we find 
Is = I c sin 4> with critical current 



I c Rd = 



64tt 




COS 



/ eV 

\2\/QvETh 
for eV > e g , T < e 

for T > e e 



(7) 



Here tt v = ttT + ^(vrT) 2 + (eV/2) 2 and 



tt/2 

i tan -1 



eV_\ 
2-nT) 



for eV > e g , T < e g 
for T > e. 



Nonequilibrium effects also influence the current- 
phase- relation Is{<p)- The rich variety of different curves 
at different voltages is displayed in Fig. 4. In contrast 
a quasi-cquilibrium theory would always predict a func- 
tional dependence as shown in the inset of the figure. 
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FIG. 4. The supercurrent-phase relation at different tem- 
peratures and control voltages. 

Since the supercurrent decays exponentially as a func- 
tion of both T and V, one might try to describe the 
system properties by a quasi-equilibrium theory with ef- 
fective I^-dependent temperature T* . Exactly this strat- 
egy was adopted in Ref. [ § with T* = y/T 2 + ~f 2 (eV) 2 . 
The the best fit to the experimental data was obtained 
for 7 ~ 6K(mV)~ 1 . On the other hand, as demonstrated 
above, the analogy between temperature and voltage is 
incomplete and even misleading. The distribution func- 
tion in the normal metal in general cannot be described 
by the Fermi function with effective temperature T* . One 
of the most striking consequences, the transition to a 
7r-junction, is not obtained in the quasi-equilibrium de- 
scription. Apart from this qualitative difference, the de- 
pendence of I C (V,T) (Q) deviates from that suggested in 
Ref. [ H . While we find a fairly good agreement between 
our results and the datafl at the temperatures of the ex- 
periment, we cannot summarize our results in general by 
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an effective temperature model. It would be desirable 
to expand the experiments to a parameter range where 
the difference between the two descriptions become more 
pronounced. 

The configuration of Fig. la is but one realization of 
a mcsoscopic system where the supercurrent can be con- 
trolled by an externally applied voltage. Another realiza- 
tion, also studied in Ref. [ || is depicted in Fig. lb. In 
this case the distribution function in the N-layer between 
two superconductors is driven out of equilibrium by the 
normal current flowing parallel to the supercurrent. Pro- 
vided d *C L there is practically no voltage drop across 
the SNS junction and only dc Josephson effect can be 
considered. Since the distribution function has the same 
form (||) as before, the previous results for Is(V) apply 
also for the structure of Fig. lb. Again good agreement 
with the experimental findingso is observed. 

Yet another system with-a. voltage controlled supercur- 
rent was studied by VolkovM. He considered a SINIS sys- 
tem where the normal metal was thin and separated by 
low transparency barriers from the superconductors. In 
this case the superconductors and the normal metal are 
in equilibrium, except that their electrostatic potentials 
are shifted relative to each other, with the total voltage 
drop across the barriers. Although this nonequilibrium 
situation is very different from the one discussed here, he 
also finds a voltage-dependent supercurrent reduction as 
well as a transition to a 7r-junction. 

If L is not large compared to ds and d, the conversion 
between super- and normal currents in the junction area 
cannot be neglected and the theoretical analysis has to be 
extended. For instance the distribution functions depend 
both Ojn_jE and y. A thorough analysis of the microscopic 
theoryQu reveals further that the expression for the cur- 
rent has to be extended, and that odd and even com- 
ponent (in energy) of the distribution function, /l and 
/t, obey two coupled diffusion equations, but with dif- 
ferent, position-dependent effective diffusion coefficients. 
They are coupled by terms of the form Inr/ e • V/l/t In 
this case the suppression of /t influences /l and weakens 
the performance of the transistor, although qualitatively 
the physical situation remains unchanged. 

We now turn to the important practical question: how 
efficient is this device as a transistor? The control and 
signal voltage, V and Vg = I c Rd, are both of the order 
of the Thouless energy E^. Thus, no voltage gain is ob- 
tained. However, the power amplification is proportional 
to the ratio of the relevant resistors Rj,jRd- Here Rr, is 
the resistance of the normal metal of length L, while Rd, 
introduced in (U), is the resistance between the super- 
conducting electrodes in a situation where I c is low and 
this transport is dissipative as well. Since both are gov- 
erned by the same material-dependent conductivity the 
ratio depends on the relevant lengths, Rl/RcL oc Ld$/d 2 . 
Hence, by choosing a sufficiently long control line, L ^> 
d,ds, a power amplification can be achieved. The limita- 
tion in operation frequency in the mesoscopic regime is 
also provided by the Thouless energy. 



In summary, we have presented a microscopic descrip- 
tion of nonequilibrium electronic properties of mesoscopic 
SNS heterostructures. The distribution function in the 
normal metal can be driven far from equilibrium by a 
voltage applied at a distance ~ L from the junction. This 
distance is limited only by the inelastic relaxation length 
l m . We analyzed how the supercurrent across the sample 
is reduced by this control voltage. The strongest reduc- 
tion and, hence, best performance of the device is found 
in a mesoscopic situation, when the distribution function 
deviates significantly from a local equilibrium form. We 
established the connection to experiments and suggested 
further tests. The possibility to control the supercurrent 
by an external voltage allows several technical applica- 
tions, for instance the use as a high-frequency transistor 
with power gain proportional to the ratio between length 
and width of the normal wire. 

We want to express our respect to the late Albert 
Schmid, who pioneered the theory on which this work 
is based. We acknowledge further discussions with R. 
Raimondi, A. Morpurgo, J. J. A. Baaselmans and H. Poth- 
ier. This work was supported by the DFG through SFB 
195 and the Graduiertenkolleg "Kollektive Phanomene 
im Festkorper", and INTAS Grant 93-790-ext. 



A. Schmid in Nonequilibrium Superconductivity, Phonons, 
and Kapitza Boundaries, ed. K.E.Gray, (Plenum, N.Y., 
1981); and further articles ibid. 

2 several articles in 'Nonequilibrium Superconductivity' eds. 
D.N. Langenberg and A.I. Larkin (Elsevier, 1986) 

3 see articles in Mesoscopic Electron Transport, eds. L.L. 
Sohn at al., (Kluwer, Dordrecht, 1997); Mesoscopic Super- 
conductivity, eds. F.W.J. Hekking et al., Physica B 203, 
201-538 (1994); C.J. Lambert and R. Raimondi, J. Phys. 
Cond. Mat. 10, 901 (1998). 

4 H. Pothier et al, Phys. Rev. Lett. 79, 3490 (1997); Z. Phys. 
B 104, 178 (1997). 

A. Morpurgo, B.J. van Wees, and T.M. Klapwijk, to appear 
in Appl. Phys. Lett. 

6 A.I. Larkin, and Yu.N. Ovchinikov, Sov. Phys. JETP 41, 
960 (1976). 

7 A.F. Volkov, A.V. Zaitsev, and T.M. Klapwijk, Physica C 
59, 21 (1993). 

8 A. A. Golubov, F.K. Wilhelm, and A.D. Zaikin, Phys. Rev. 
B 55, 1123 (1997) and references therein. 

9 F.K. Wilhelm, A.D. Zaikin, and G. Schon, J. Low Temp. 
Phys. 106, 305 (1997). 

A.A. Golubov and MYu. Kupriyanov, Zh. Eksp. Teor. Fiz. 
96, 1420 (1989) [JETP 69, 805 (1989)]; W. Belzig, C. 
Bruder, and G. Schon, Phys. Rev. B 54, 9443 (1996). 

1 F. Zhou et at, J. Low Temp. Phys. 110, 841 (1998). 

2 L.N. Bulaevskii, V.V. Kuzii, and A.A. Sobyanin, Solid St. 
Comm. 25, 1053 (1977). 

3 H.Courtois et al, Phys. Rev. Lett. 76, 130 (1996). 

4 A.F. Volkov, Phys. Rev. Lett. 74, 4730 (1995). 



4 



